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Vorticity aligned cylindrical floes of carbon black 
particles are formed in steady flow at low shear rates, 
and strikingly, appear as transient structures in the flow 
response of gels produced by the quenching of high rate 
shear thickening flows. 

Structure formation in response to flow is common in 
soft materials and complex fluids and is dependent on 
the nature of the cohesive forces in the system, time 
scales for diffusive motions and hydrodynamic interac- 
tions. For Brownian systems with hard sphere or repul- 
sive interactions, distortion of the microstructure due to 
an imposed shear rate occurs when the timescale for flow 
is smaller than that for diffusion. The relative strength 
of flow forces compared to thermal or diffusive forces can 
be gauged by the dimensionless Peclet number 



For Pe <C 1 the system displays isotropic arrangements 
of particles, but for Pe > 1 the timescale associated 
with flow dominates that for diffusive relaxation. Under 
these conditions, the system becomes anisotropic due to 
separation of particles along the elongational axis and 
aggregation along the compressional axis [1] . String for- 
mation by alignment along the flow direction is typical 
in dilute systems at high shear rates and further struc- 
turing may take place by periodic arrangement of these 
strings in the vorticity-gradient plane. In more concen- 
trated suspensions, hydrodynamic interactions lead to 
the formation of hydro-clusters and a shear thickening or 
jamming transition [21 Ej. In systems with constitutive 
instabilities due to a negative slope in the stress-strain 
curve, co-existence between different deformation states 
under a single applied stress gives rise to shear banding 
along the velocity gradient direction [3]. For attractive 



systems, structure formation as discussed above is also 
found. Additionally, the alignment of particles or con- 
stituents along the vorticity axis has been reported and 
has emerged as a somewhat general phenomenon [5j |6] . 
Shear thickening though, is not generally observed in 
colloidal systems with sufficiently strong attractive in- 
teractions to form flocculated gels and is presumed not 
to occur [7J . It has been reported in a handful of systems 
to date [HI O HOI , although the complex composition of 
the suspensions studied in some cases makes it difficult 
to unambiguously study the mechanism. 

We have studied the steady state flow behavior of dilute, 
simple hydrocarbon dispersions of carbon black particles 
and, surprisingly, observe shear thickening above a com- 
position dependent critical flow rate, % « 10 2 — 10 3 s -1 . 
The shear modulus of gels formed by pre-shearing above 
the critical shear thickening rate displays an interesting 
power law dependence on the stress applied during the 
pre-shear. This is well accounted for by the dependence 
of the cluster size on the applied shear stress, and the 
resulting increase in cluster number density [IT] . Subse- 
quent deformation of these shear thickened gels at low 
shear rates produces highly oriented vorticity aligned 
floes that gradually break down into small isotropic clus- 
ters over timescales of « 300 seconds. 

We use tetradecane (Aldrich Chemical Co.) dispersions 
of 2 to 8 wt.% of 0.5 fim carbon black particles (Cabot 
Vulcan XC72R) with a fractal dimension df = 2.2. Rhe- 
ological data reported here were recorded using an AR- 
G2 (TA Instruments) rheometer in strain control, al- 
though separate experiments were also run for verifica- 
tion purposes on a pure- strain controlled ARES-LSI in- 
strument. A variety of geometries, standard and rough- 
ened, were used to check for and to mitigate the impact 
of possible wall slip. Although from our measurements 
it appears that wall slip is not significant, its preva- 



lence in other studies on "sticky" colloidal gels makes 
the interpretation of data at low shear rates trouble- 
some [12]. A more detailed study of the low shear rate 
regime is underway. Optical observations of microstruc- 
tures under shear were made using a Bohlin CS rheome- 
ter with a transparent base-plate through which images 
were recorded (by reflection off the rheometer top tool) 
using a CCD camera. 

The flow curves for samples in steady shear display 
thixotropy, as is typical for flocculating systems. Steady 
state flow curves of the shear rate dependent stress dis- 
play a composition-dependent, continuous shear thicken- 
ing transition at high shear rates, 7 C « 10 2 — 10 3 , with 
distinct increases in the stress and thus the viscosity 
of the dispersions with increasing shear rate, Figure [2j 
From optical observations, Figure [TJ at low shear rates 
7 < Is -1 , the system is composed of large pieces of the 
fractured gel which gradually break into smaller pieces 
with increasing shear rate. Around 7 ~ 10° — lO 1 ^ -1 , we 
observe a slight, if irregular increase in viscosity with in- 
creasing shear rate. Coincident with flow in this region, 
we observe the aggregation of clusters along the vor- 
ticity axis, and the formation of elongated floes which 
roll with the shear flow. Such floes are also observed in 
long observations overnight at constant shear rate, in- 
dicating that, within the experimental time frame, they 
represent the steady state response of the system to flow 
in this range of shear rates. Continued increase of the 
shear rate results eventually in dissolution of the cylin- 
drical structures into isotropic clusters. These clusters 
densify until 7 C beyond which there is a transition to 
the shear thickening flow as breakage of the aggregates 
is observed, consistent with an increase in the effective 
volume fraction of particles in the dispersion. 

Previous studies have implicated negative first normal 
stress differences, Ni = an — 022, in the formation of 
vorticity aligned structures in flocculating systems [SUB]. 
Normal stresses arise due to elasticity in the flowing sus- 
pension with elongational stiffness and bending rigidity 
both contributing p~3l [111 [15] • I n our system, we are un- 
able to identify any systematic changes in normal stress 
at low shear rates. It may simply be that the stresses 
are too small to be accurately measured in our rheome- 
ter. In the shear thickening regime however, samples 
develop large negative normal stresses, on the order of 
Ni « 2— 10.(T s h ear , depending on composition, due to the 
large inertia developed at these high shear rates. Abrupt 
cessation of shear in this regime produces a quench into a 
shear-thickened gel state where, remarkably, both shear 
stresses and normal stresses were found to persist for 
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Figure 1: Microstructure under shear in parallel plate 
geometry with a gap size of 100 jim for 3 wt.% sample, 
(a) 7 = 6.67s- 1 40 mm field of view showing cylindrical 
floes aligned along the vorticity direction, (b) 1.5 mm 
field of view, 7 = 133s -1 during shear (c) 1.5 mm field 
of view, 7 = 1330s -1 during shear. To obtain sufficient 
light transmission through the optically dense sample in 
(c), the illumination was set several times higher than 
that used for (b). Binning of pixels on the CCD was 
used to decrease the required exposure time and so the 
resolution of image (c) is half that of (b) . 
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Figure 2: Steady state flow curve of carbon black disper- 
sions from 2 to 8 wt. %.+ : 2%, > : 3%, o : 4%, o : 5%, 
A : 6%, □ : 8% The dotted-line shows the stress due to 
the background viscosity of the solvent. 



long times (> 1 hour) under quiescent conditions and 
to dissipate rapidly on application of shear flow. We 
have studied the dynamics of the relaxation of the inter- 
nal shear stresses and find that they relax with a power 
law dependence U{ ~ t~ 0A over timescales as long as 
10 3 — 10 4 s [TT] . Accurate measurement of the modest 
normal stresses in dilute suspension rheology is made 
difficult by the lack of better than « 10-20 Pa resolution 
available on most commercial rheometers, which rely on 
a strain gauge style mechanism that is usually not stabi- 
lized against temperature fluctuation. The presence of a 
large inertial contribution N[ nertial = -3/20 puj 2 R 2 
to our signal further complicates this in our case, but 
from more concentrated samples, (j) > 4 wt.%, we could 
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reliably extract the sample contribution to the normal 
stress, Figure [3} Correction of the data for inertial con- 
tributions shows that a positive normal stress difference 
is developed at the onset of shear thickening, consistent 
with the idea that shear thickening, as indicated from 
optical observation, is due to cluster breakage and an 
increase in the effective volume fraction of the system. 
The negative normal stress observed at rest relative to 
the flowing state appears to originate during the rapid 
quench from the tendency of the particles to collapse 
into dense clusters on cessation of shear. Accurate mea- 
surements of the normal stress transients as well as the 
normal stress dynamics in the quiescent state are the 
subjects of continuing work. 
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Figure 3: First normal stress difference (relative to start- 
ing value, 7Vi(0))for various compositions as a function 
of shear rate. Data are averaged over 4 runs. Error bars 
correspond to one standard deviation. 



The deformation of quenched shear thickened gels by 
steady flow (0.1s -1 < 7 < 10s -1 ) or large ampli- 
tude oscillatory shear (u = 1 rad/s, 7 ^ 50%) results 
in the rapid emergence of highly anisotropic vorticity 
aligned structures. The transition from shear thickened 
3-dimensional gel to an assembly of discrete cylindri- 
cal floes in rolling flow for a sample of 2 wt.% carbon 
black, using a parallel plate geometry, is shown in Fig- 
ure [4] (see supplemental data for video). In contrast to 
the structures observed in the steady state response at 
low shear rates, these transient floes were much more 
sharply defined and with aspect ratios exceeding 10 2 , 
spanned macroscopic portions of the shear cell. 
Fourier transformations of the images provide a means 
of following the formation and alignment of the vorticity 
aligned floes. We inspect a small section of the sample, 3 
mm on a side and 2.5 mm from the tool edge, to limit the 





Figure 4: Carbon black dispersions, (p = 2 wt.%, shear 
thickened by strong flow in the secondary thickening 
regime. Samples are quenched to zero-shear rate and 
then subsequently deformed, starting at t = with a 
shear rate of 7 = 10s -1 . Parallel-plate geometry with 
gap-size d = 50/xm. Field of view is 1.35 cm in width. 
Vorticity direction is along the bottom-left to top-right 
diagonal as given by the arrow in the first image. 



distortion due to the circular geometry. We are unable 
to resolve the rapid changes in the characteristic length 
scale of the system due to the low optical resolution of 
the camera as well as the short timescale at which a 
stable structure emerges, Figure [5| We can, however, 
monitor the development of alignment in the system. We 
extract the azimuthal intensity dependence of the FFTs 
at the spatial frequency corresponding to the periodicity 
of the floes and model them with a Gaussian. As shown 
in Figure [5j alignment is rapid and peaks around 20s, 
7 = 200, after which it quickly decreases as the floes 
break up. 

Using the parallel-plate geometry, we examined the de- 
pendence of the floe geometry on physical confinement. 
The width of the rolling log-like structures measured at 
a distance of 2 mm from the tool edge, wi og , was moder- 
ately larger than the gap size, d, and varied in direct 
proportion with the gap, Figure [6j The accuracy in 
determining the log width is limited by the 20/xra per 
pixel resolution at which data was taken, but in general, 
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Figure 5: (a) 3 mm sections of optical micrographs over 
time. Gap = 50fim. From left to right: t — 1,3, 12,24 
and 55s (b) Azimuthal spread of intensity from FFTs of 
images in (a). Solid lines are Gaussian fits to the data. 
Inset: Peak widths vs. log time 



wi g ~ lAd. At gaps larger than 750 fim. the structures 
could no longer be clearly discerned. From the images, 
it can be seen that smaller gaps promote sharper inter- 
faces between the structures and the surrounding fluid. 
The width of the structures increases slightly with de- 
creasing shear rate in moving from the outside of the 
parallel-plates towards the center. Increasing gap size 
results in greater areal coverage of the parallel-plates as 
the periodicity of the logs does not increase commensu- 
rate with their width. This implies that the floes become 
less dense internally with increasing thickness. At long 
times, the dissolution of the structure proceeds from the 
outside edge of the tool, where the shear rate is highest. 
In observations employing a low angle cone (1°) where 
the shear rate is now constant across the tool, dissolu- 
tion proceeds from the outside edge where the gap is 
largest, underscoring the importance of confinement in 
stabilizing the floes. 

The frequency and strain dependence of the elastic mod- 
ulus of the log-like floes differ from those of the shear 
thickened gel. Here, the vorticity- aligned structures 
were produced by stopping shear flow at t = 20s in a 
40 mm, 1° cone-plate cell with an edge gap of 350 \±m. 
Both the gel and the aligned structures show strong elas- 
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Figure 6: Microstructure of 2 wt. % carbon black disper- 
sions at « 25s, 7 = 10s -1 after shear thickening flow in 
40mm plate-plate geometry with various gaps, d. Field 
of view is 1.35cm in width, (a) d = 25/xm (b) d = 50fim 
(c) d = 100/xm (d) d = 250/xra 



tic responses, with nearly frequency-independent mod- 
uli. The yield strain, however, is markedly lower for the 
structures than for the gel, Figure [T^i. Yielding of the 
vorticity- aligned structures proceeds without the pro- 
nounced upturn in the loss modulus, G" that typically 
signifies yielding in the gel, and the crossover between 
G' and G" occurs at higher strains than is typical for 
the gel. The existence of a yield strain implies that dis- 
placement and deformation in the system is not accom- 
modated solely by the rolling motion of the floes. 

Vorticity aligned structures have been observed in a 
variety of complex fluids such as thixotropic clay gels 
[IT], nanotube suspensions [6J and attractive emulsion 
droplets [5], associated in each case with negative normal 
stress differences. In our system, the carbon black dis- 
persions appear to display small negative normal stresses 
over all shear rates before increasing at the shear thick- 
ening transition and decreasing again beyond the peak 
of the transition. No distinct transition in normal stress 
was observed under steady flow at low shear rates, 
though this is subject to the difficulty involved in ac- 
curate measurements of small changes in these stresses 
due to instrument resolution. Measuring the absolute 
value of the normal stress is further complicated by the 
suction of the tool surface by the sample, i.e. in such 
systems it is not clear what constitutes absolute zero. 
Rather, significance can only be ascribed to differences 
in the normal force signal with respect to some initial 
condition. 

The rapid formation and the high degree of coherence 
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Figure 7: Strain dependent moduli of shear thickened 
gels (open symbols) and log-like floes (filled symbols) for 
carbon black dispersions of various compositions. Strain 
data taken at u = 1 rad/s. and frequency data at 7 = 
0.05%. o: 2 wt.%; A: 3 wt.%; <: 4 wt.%.(a) Storage 
modulus. Inset: Frequency dependent data, (b) Loss 
modulus. 



of the structures formed in the transient flow response 
of the gel suggest strongly that relaxation of quenched 
normal stresses play a significant role. The dissipation 
of the normal stress with continued flow precludes long 
term stability of the structures and as a result they exist 
only transiently. Our results are similar to observations 
of an elastic instability associated with flow-induced 
clustering in semi-dilute non-Brownian nanotubes [T8] . 
Here, the transient rheological response of a flowing 
suspension of nanotubes subjected to a shear-rate 
quench is underlined by the formation of vorticity- 
aligned aggregates. The existence of these aggregates 
can be parameterized by confinement and shear rate. In 
our system, however, we also observe these structures 
as a response to steady flow at low shear rates. By 
stopping the flow of shear-thickened systems during 
the transient response, we have been able to access 
measurements of the frequency and strain dependence 
of the shear modulus of these structures. The transition 
from a three-dimensionally connected network to a 
two-dimensional system of highly anisotropic structures 
and the dissolution of these structures thereafter shows 



interesting features. Locally, one can discern changes 
in topology as the floes make and break connections 
with neighboring structures, and track changes in the 
order parameter as shear finally dissolves the system 
into non-aligned clusters. The instrument used in this 
study did not permit the consistent application of 
arbitrarily small deformations. As a result, we were 
unable to conduct measurements at the low strain rates 
that would have permitted better resolution of the 
dynamics of structure formation and dissolution in our 
samples. Future work will benefit from a more capable 
instrument that also permits measurements of normal 
forces and will pursue more detailed studies of the effect 
of confinement and composition on the mechanics of 
these structures. 

The authors gratefully acknowledge C. Kim and 
H. Wyss for very helpful discussions, and the Inflneum 
Corporation for funding. 
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